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The process of baryogenesis through the evaporation of black holes formed at the end of inflation 
phase is considered. The increase of black hole mass due to accretion from the surrounding 
radiation after the reheating is taken into account. It is shown that the influence of the accretion 
on the baryogenesis is important only in the case when the initial values of black hole mass 
are larger than ~ 10'' g. The behavior of calculated baryon asymmetry, as a function of model 
parameters, is studied. 
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I. INTRODUCTION 

First discussions about possible connection between 
baryon asymmetry and primordial black hole evapora- 
tions appeared at the middle of seventies, just after dis- 
covery of a phenomenon of the black hole evaporation. 
The possibility of an appearing of the excess of baryons 
over antibaryons in the process of evaporation of pri- 
mordial black holes was noticed in works The 
mechanism discussed in |^J^ doesn't require the non- 
conservation of baryon number in an underlying micro- 
scopic theory. Detailed calculations using this mecha- 
nism were done in works . 

Later, this question was studied in the context 
of grand unified theories (GUTs) (i.e., theories in which 
baryon number is not conserved). The idea is simple: 
when black holes decay by the emission of Hawking radi- 
ation, they may emit baryon-number- violating Higgs par- 
ticles (and/or leptoquarks) whose decays naturally gen- 
erate baryon asymmetry. 

In the work of J.Barrow et al. |^ baryogenesis via pri- 
mordial black holes was considered using GUT and ex- 
tended inflation scenario. The formation of very light 
primordial black holes (which disappear, as a result of 
evaporation, before nucleosynthesis, without any trace, 
except of the net baryon asymmetry) seemed to be most 
probable in inflationary models with first order phase 
transition. Recently, however, it was shown that in some 
modern variants of second order inflation models |^,^ the 
formation of small black holes (and, consequently, baryon 
asymmetry production via black hole evaporation) is also 
quite possible. 

In two recent works the baryogenesis through the 
evaporation of primordial black holes was studied quan- 
titatively. Authors of these works argued that this sce- 
nario of baryogenesis can, in principle, explain the ob- 
served baryon number of the Universe (which is con- 
strained by primordial nucleosynthesis data to be in 
range (1.55 — 8.1) -10^^^). However, the distinct conclu- 
sion of works [^,0 is that a sufficient number of black 
holes can survive beyond the electroweak phase transition 
and therefore, the baryon excess produced is not washed 
out by sphaleron transitions. These conclusions are in 
contradiction with qualitative statements of pioneering 
work of J.Barrow et al. Authors of ||^,|l^ claim that 



the difference between their result and those of J.Barrow 
et al. can be explained by taking into account the ac- 
cretion (accretion term was omitted in formulas of work 
§)• 

Aside from the problem of the possible erasing of pro- 
duced baryon asymmetry by sphaleron transitions, it is 
important to have the definite answer on the following 
question: is it possible to obtain by such mechanism (i.e., 
by primordial black hole evaporations) the cosmologically 
interesting value of baryon asymmetry? 

In present work we try to answer just this question. 
We calculate the baryon asymmetry using the same as- 
sumptions, as in the work of J.Barrow et al. The main 
attention is paid to rigorous solution of kinetic equations 
(containing accretion term) and to comparison of exact 
results with predictions based on approximate formulas 
derived in fql. 



II. BASIC ASSUMPTIONS 

1. A generation of the black holes in early Universe 
may take place in many inflationary scenarios. The main 
assumption is that in the early Universe there was the 
period of inflationary expansion and the inflation is com- 
pleted by a symmetry-breaking phase transition. The en- 
ergy scale of symmetry breaking is co ^ 10^^ GeV (GUT 
scale). 

Clearly, the concrete mechanisms of black hole produc- 
tion are not the same in different inflationary models. 

la. Inflation of "old" type (extended inflation, flrst 
order inflation): the Universe exits from a false- vacuum 
state by bubble nucleation. Reheating and thermaliza- 
tion of the Universe proceeds through bubble collisions 
at the end of inflation ||^ . False vacuum energy density 
is 



so 



py - lO'"' GeV^ = Ar. 
Hubble parameter at the end of inflation is 
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so the time of the end of inflation is 
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Usually one supposes that 



end 
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tf < Th 



(2.4) 
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{tRH is the time of reheating, i/ is the time of formation 
of black holes, is life-time of black hole). 

Black hole production proceeds 1) via the gravitational 
instability of inhomogeneities formed during the thermal- 
ization phase (i.e., during the bubble wall collisions) or 
2) via the appearance of trapped regions of false vacuum 
caught between bubbles of true vacuum. 

lb. "New" or "slow-roll-over" inflation: in this case 
black holes are produced through collapses of the over- 
dense regions in space. So, for large probability of pri- 
mordial black hole production one must exist large ampli- 
tudes of primordial density fluctuations at small scales. 

Such fluctuations appear during inflation. The over- 
dense region of mass M can produce black hole when 
this fluctuation crosses horizon inside. At this time M is 
equal to horizon mass M/i, and black hole produced has 
mass Mbh which is close to M/j, 



Mbh - O.IM,, 



(2.6) 



The time of PBH formation is 
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rripi nipi 
If, e.g., Mbh is about 10^ g, then 
tf - 10'^^ s. 
At this time, the corresponding scale factor is 

a/ ~ (i/o^'(2.4 • 104)-3/4t3/2)l/3 _ 10-27 

and the comoving length scale of the perturbation is 

Ctf 



(2.7) 
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10-14 PC. 



Further, we use the known formula 
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3 lOi^GeV 



(2.10) 



1, Trh 
3 lOiOGeV 



(2.11) 



Here, M ~ 10^^ GeV (it is ^/p^), Trh is reheating tem- 
perature; Nf is a number of e-folds before the end of in- 
flation beginning from the moment when the scale crosses 
horizon outside. 

If, e.g., Trh ~ 10" GeV, one has N),^ ~ 1. Now it is 
clear that the fluctuations responsible for Mbh ^ 10^ g 
are formed just near the end of inflation. So, in this 
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FIG. 1. The time evolution of scale factor parameter a, for 
/3 = 0.1 and for two values of Mq. 
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10"^ g are lightest 



sense, black holes with Mbh 
ones. 

The large amplitudes of density perturbations at small 
scales corresponding N\j 1 are naturally obtained in 
the hybrid inflation model of refs. 0,^. 

If TiiH ^ 10^"'^ GeV, the reheating time is 



tRH = 0.3- 
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(2.12) 



So, black holes can be produced even at a time before 
reheating, when the Universe is dominated by the oscil- 
lations of the inflation field. 

2. We assume that all produced primordial black holes 
have the same mass. It will be approximately so in hy- 
brid inflation-type models, where the sharp maximum 
of density fluctuation amplitude exists at some definite 
scale. 

3. We assume that at tjin the part of energy density 
is in black holes: 

p{tRH) - l3pBH{tRH) + (1 - l3)pR{tRH)- (2.13) 

j3 is the free parameter of our model. 

In the following we will consider only black holes hav- 
ing the life-time Th, which is much larger than t^H, 



Th > tRH- 



(2.14) 



We will see that the most interesting predictions for 
baryon asymmetry don't depend on tjin and are the same 
for both inflation scenarios provided we use in both cases 
the parameter [3. 
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FIG. 2. The dependence of function R on r, for /3 = 0.1 FIG. 3. The r - dependence of the temperature of radiation 

and for two values of Mo. in the Universe, for /? = 0.1 and for several values of Mq. 



III. KINETIC EQUATIONS 

Evolution of black hole mass is described by the equa- 
tion 



(3.1) 



Accretion term AMgj^ is equal to crabs Prc, where aats 
is cross section of absorbtion of relativistic particles by a 
black hole, 



q2 

(Jabs = 27T:Mgu ■ — , 



(3.2) 



pii is an energy density of the radiation. Evaporation 
term is — where a (to) counts the degrees of freedom 
of the black hole radiation (m is an instantaneous value 
of the black hole mass) . 

In the following we will use the value P^,E3| 



a{m) = const = 80 • 10^^ — . 



(3.3) 



Evolution of the radiation energy density is given by 
the equation 



PR 



-4:-pR - MBHriBH, 



nBH{t) = 



PBH^tRH) a^{tRH) 



Mo 



PBH{tRH) = Pp{tRH)- 



(3.4) 

(3.5) 
(3.6) 



Here, ubh is number density of primordial black holes, 
Mo is initial value of black hole mass, a{t) is scale factor. 

Evolution of scale factor is given by Friedmann - Ein- 
stein equation: 



SttG 



{Pi 



M 



BH 



riBH)- 



(3.7) 



In this paper we use, for simplicity, the Einstein - de Sit- 
ter model. 

Equation for evolution of the baryon number is 

Mb//c2 a 

'^B = -iTT7^ nBH ■ efn - S-UB- (3.8) 

6k1 bh o, 

We assume here, approximately, that for cLAIbh change 
of black hole mass we have (IMbh ■ c? /ikTBH particles 
produced {Tbh is temperature of black hole connected 

with Mbh by Hawking relation kTBH = sttGMbh 
is the fraction of X particles decaying with violation of B. 
Typically, fn ^ ^ ^ 10~^. Baryon number is obtained 
by introducing the factor e which is given by the relation 

,^y2B.m^Ill_Iil^. (3.9) 
I X 

e 7^ if C,CP are not conserved. For obtaining baryon 
asymmetry one must divide ub on entropy s which is 
given by 



Svr^ ^, , f k 
.W = ^...T3(0 - 



(3.10) 



It is very convenient nM to use dimensionless variables: 
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FIG. 4. Baryon asymmetry parameters as a function of Mo 
for three values of l3 (solid lines). Dashed lines: calculation 
by approximate formula ([l.2|). Dotted lines: calculation by 
the formula (4.2) with substitution (i.ll). 
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(3.12) 



for which we have simple initial conditions 

■m{T = tjih) = a{T = trh) = R{t = trh) = 1, 



(3.13) 



TRH — V Gp{tRH)tRH- 



We assume, approximately, that X-particles are evapo- 
rated when the following condition holds: 



SkTsH = 3/c- 



> Mxc^ 
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BH 



SttMx ' 



Here, Mx is a mass of X particles. 

Finally, we have the following free parameters: 
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f3 
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(3.14) 



(3.15) 



(3.16) 



Some re sults o f th e s olut ion of the system of kinetic 
equations (3J), (3^), ( |3.7D and ( |3.8| ) are presented on 
Fig. 1-4. All calculations are carried out with pitun) ~ 
10^^ GeV''. Such value of p{tRH) corresponds to Irh ^ 
10"'^^ s and reheating temperature Trh ~ 10^'^ GeV. 



IV. APPROXIMATE FORMULAS FOR B 

One can easily show |^ that for practically important 
case, when 



/3> 



tRH 
Th 



(4.1) 



and when accretion is not important (Mq < lO'^g), there 
is an approximate solution of the kinetic equations for B 
which is given by the formula 
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V^(-] p"HtRH)y. 



m^iMo g* \ Th J 



(4.2) 
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Using the relation 



m 



p{tRH) - ■ -p 72—' 
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(4.3) 



one can see that B doesn't depend on pitun) or t^H- 
The expression for r is obtained from the equation 



Mbh = - 



and is given by 



Th = 



Ml 
3a 



For Mq ~ 10^ — lO-^ g one has 

TH ~ 10-21 _ 10-18 s. 

In our numerical calculations we used t j^jj 
one has 

tRH < T/i • 



(4.4) 

(4.5) 

(4.6) 
10-33 s. So 



(4.7) 



From here, one has the condition for parameter (3: 



/3>^/— -(lO-^-lO-^- (4.8) 

For such values of /3, tun, Th one has the final formula: 

f Mo Y'^' 

■J ^ 

(4.9) 



B ^ IQ- 



Vio3 

Mx 
10i4GcV 



Condition (4.1) means M that the black hole energy 



density dominates at the time of evaporation. If the value 
of parameter f3 is such that 



/3< 



tRH 
Th 



(4.10) 



the evaporation occurs while background radiation dom- 
inates the energy density of the Universe (with the en- 
tropy arising mainly from the radiation). In this case 
the approximate formula for B ^ can be obtained from 
Eq.(4.2) by substitution of last factor: 



\ThJ ^(l-/?)3/4- 



(4.11) 



i s we ll described by the ap proxi mate formulas (4.2) and 
(4.2) with substitution of ( 4.11 ). These formulas were 
obtained without taking into account the accretion, and 
it means that the accretion process is not important for 
small values of initial black hole masses. At larger masses 
accretion becomes to be important (the same was argued 
in the work ||l^). One can see from Fig. 4 that the accre- 
tion leads to a significant decrease of B (this decrease is 
a consequence of a significant growth of entropy after an 
evaporation of black holes enlarged due to accretion). 

The resulting formula for B (Eq.l.S) doesn't depend 
on pitfin) and on tun and, therefore, is valid for both 
variants of inflationa ry s cenario mentioned in Sec.|| (as 
far as the condition ( [4.lD is fulfilled). It is seen that, if 
Mo ^ 103 g and Mx ~ lOi^'GeV, the predicted asym- 
metry is quite small (even for e ~ 1). This our calcu- 
lation strongly disagrees with the corresponding results 
of work | [Tc[ |. The effect (baryon asymmetry) can be 
large if Mq < 10^ g and/or Mx < 10^'' GeV (e.g., if 
Mx ~ 10" GeV). 

It follows from Fig. 3 that the temperature of the Uni- 
verse at a moment of the evaporation is smaller than 
^ 100 GeV only in the case, when Mq > 10^ g. For 
smaller values of Mo (for which the value of B can 
be acceptably large) one has T^v > 100 GeV, and the 
problem connected with the sphaleron transitions exists. 
Evidently, the baryon asymmetry produced by primor- 
dial black hole evaporations can survive only if nonzero 
{B ~ L)-value is generated in decays of Higgs particles of 
GUT. 
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V. CONCLUSION 



The main conclusion of our work is the following: the 
predicted baryon asymmetry in the region Mq < 10^ g 
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